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A recently developed linear-scaling density- functional theory (LS-DFT) formalism is used to calcu- 
late optical absorption spectra of hybrids of Ceo and the conjugated polymers poly(para-phenylene) 
(PPP) and poly(para-phenylene vinylene) (PPV). The use of a LS formalism allows calculations 
on large systems with realistic proportions of Ceo, which has been of interest for the use of such 
materials in photovoltaics. Two different bonding structures are tested for the hybrid PPP and 
for both systems additional peaks are present in the absorption spectra below the original onset of 
absorption. By identifying the eigenstates involved in the relevant transitions, a weighted density 
difference is formed, demonstrating the transfer of charge between the polymer chain and the Ceo, in 
agreement with experiment. For the hybrid PPV, no additional peaks are observed in the absorption 
spectrum. 



I. INTRODUCTION 

Since the discovery in 1976 that conjugated polymers 
could be made conducting^—, they have been the subject 
of much research, both for reasons of fundamental theo- 
retical interest and due to the promise of a wide range of 
exciting applications. They are highly attractive candi- 
dates for use in applications ranging from light-emitting 
diodes (LEDs) to photovoltaic cells^— . This is in part 
due to the ability to influence their properties via doping 
and modification of side chains, enabling for example the 
creation of materials with a wide range of energy gaps 
resulting in materials ranging from metallic to insulat- 
ing 2 -^. Furthermore, they benefit from useful mechanical 
properties such as the ability to be flexible, and from 
desirable processing methods, particularly the ability to 
process from solution, which also has the important ef- 
fect of reducing the final cost of the materials^™. As a 
result, a great deal of effort has been extended into de- 
veloping and studying such polymers with the hope of 
replacing inorganic semiconductors in some existing ap- 
plications as well as creating new ones. 

For the purposes of this work, we are particularly in- 
terested in the application to photovoltaic cells, where 
polymer based devices could have a number of advan- 
tages over the standard silicon based solar cell, particu- 
larly in terms of cost. Indeed, conjugated polymers such 
as poly(para-phenylene) (PPP) and poiy(para-phenylene 
vinylene) (PPV) (whose structures are depicted in Fig.QJ 
have been suggested as possible candidates for use as 
a donor when in combination with an acceptor such 
as the fullerene Cgo^- For this use, certain criteria 
must be met, including the occurrence of charge transfer, 
which has been demonstrated in fullerene-polymer mix- 
tures^—. Indeed, Ceo-PPV photovoltaic cells have been 
realized in a layered structur o 13 i 14 . 

A useful way of studying these systems and compar- 
ing the behaviour of different conjugated polymers is 




(a)poly(para-phcnylcne) 




(b)poly(para-phenylene vinylene) 

FIG. 1. Schematics showing the structure of segments of PPP 
and PPV, with C atoms shown in white and H atoms in black. 



through computational simulations, in particular first- 
principles calculations which allow the calculation of op- 
tical absorption spectra, among other useful properties. 
One very popular methodology is density-functional the- 
ory (DFT ) 15 i 16 , which has been applied extensively to a 
large number of systems and shown to give remarkably 
accurate results for ground state properties, and can also 
be used as a good first approximation to excited state 
properties such as absorption spectra. The accuracy of 
these spectra will of course be limited by the approxi- 
mations inherent in practical implementations of DFT, 
and arising from the fact that the Kohn-Sham eigenval- 
ues have no rigorous correspondence with the true quasi- 
particlc energies. Indeed, one can only achieve good 
agreement with experiment by applying the "scissor op- 
erator", where the conduction states are rigidly shifted 
upwards in energy to correct for the underestimation of 
the band gap by DFT. However, previous calculations on 
PPP have shown good agreement with the experimental 
absorption spectra using the local density approximation 
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(LDA)** and so we expect the use of DFT to be a rea- 
sonable approximation for these systems. Furthermore, 
the systems of interest here involve many atoms, even 
when limiting ourselves to e.g. a single chain combined 
with a single molecule of Ceo, so that the cost of more 
accurate methods such as the GW approximations- 
is prohibitively expensive for calculations on systems of 
this size. Indeed, even conventional DFT codes can be- 
come too expensive due to the cubic scaling bottleneck 
encountered for systems over a given size, and so it be- 
comes useful to have a linear-scaling (LS) formulation 
of DFT, which takes advantage of the principle of near- 
sightedness and therefore allows one to access truly large 
systems. 

To this end we have employed a recently developed 
methodology for calculating first the Kohn-Sham con- 
duction states and then the imaginary component of the 
dielectric function in a LS-DFT framewor k 21 ' 22 . This 
new methodology combines the optimization of a sec- 
ond set of localized orbitals to describe the conductions 
states with the use of Fermi's golden rule and has been 
implemented in the LS-DFT code oneteb* 3 .*^*. We have 
applied this methodology to the study of Cgo hybrids of 
PPP and PPV. 

In this paper we first outline the basic theory behind 
the ONETEP code, including the method used to calculate 
the Kohn-Sham conduction states. We then describe the 
method applied to the calculation of optical absorption 
spectra. Finally we will present results for the absorption 
spectra of the hybrid PPP and PPV, before summarizing 
our conclusions. 



II. METHODOLOGY 

ONETEP is a linear-scaling DFT code which takes ad- 
vantage of the principle of near-sightedness, and in par- 
ticular the exponential decay of the density matrix (DM) 
in systems with a band gap**~— . This is achieved using 
a DM formulation, within a strictly localized represen- 
tation of non-orthogonal generalized Wannicr functions 
(NGWFs)**. The NGWFs are themselves represented in 
terms of a periodic cardinal sine (psinc) basis set**, which 
allows them to be optimized during the calculation. The 
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DM is written in separable for m 3 as: 

<r|p|r')=X> Q (r)A-^ (r') 

a/3 



(1) 



where K a " is the so-called density kernel and {4> a (r)} 
are the NGWFs. An energy cut-off is applied to the 
psinc basis set in an analagous manner to plane- wave ba- 
sis sets**, and thus ONETEP combines the advantage of a 
systematically improveable plane- wave basis set with the 
benefits of a localized basis set. In a ground state energy 
calculation, both the DM and NGWFs arc optimized to 
convergence, with the density and potential updated un- 
til self-consistency is achieved in the usual manner. 



During the calculation the eigenstates are not directly 
referenced, however they can be recovered via a single 
diagonalization at the end of the calculation. This is an 
0(A?" 3 ) process, however the prefactor is small as we op- 
timize a minimal set of NGWFs that is much smaller 
than the number of psincs, so that the cost is negligible 
compared to the overall simulation time. This enables 
for example the calculation of densities of states (DOS) 
and band structures. However, as the NGWFs are op- 
timized to minimize the total ground state energy, they 
are not expected to accurately represent the Kohn-Sham 
conduction states, and in some cases certain conduction 
states are found to be completely absen t 21 ' 22 ' 34 . In or- 
der to be able to calculate optical absorption spectra one 
therefore needs a method of optimizing a new set of NG- 
WFs to accurately represent the unoccupied states. This 
is achieved by defining a projected Hamiltonian: 
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where p is the density operator and a is an energy shift. 
The lowest eigenstates of this projected Hamiltonian are 
the lowest conduction states of interest and so it can 
be used in a similar procedure to ground state ONETEP 
calculations, optimizing a second set of NGWFs for a 
fixed number of conduction states, omitting the self- 
consistency process and instead re- using the ground state 
density and potential. Following the conduction NGWF 
optimization, the valence and conduction NGWFs are 
combined to form a joint valence-conduction basis set, 
which is capable of representing both the occupied and 
unoccupied Kohn-Sham states. Further details have been 
given elsewher o 21 ' 22 , including a demonstration of the 
linear-scaling behaviour of the method. 

Once one has a process for calculating the Kohn-Sham 
conduction states as well as valence states, Fermi's golden 
rule can be applied to the calculation of optical absorp- 
tion spectra. The methodology** used here is the same 
as that applied in the cubic-scaling plane- wave pseudopo- 
tential (PWPP) code CASTEP**. Using the dipole ap- 
proximation, whereby only first order terms are included, 
one can write the imaginary component of the dielectric 
function as: 
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(3) 

where v and c denote valence and conduction bands re- 
spectively, IV'k} is the nth eigenstate at a given k-point 
with a corresponding energy E£, f2 is the cell volume, q 
is the direction of polarization of the photon and hio its 
energy. This expression can be directly applied for molec- 
ular systems but for extended systems one must use the 
momentum rather than position operator, as the position 
operator is ill-defined for periodic boundary conditions. 
The correct relation between the two is**: 



V nh r 



(4) 
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where it is particularly important to include the second 
term, which involves the commutator between the non- 
local pscudopotcntial and the position operator. In this 
manner one can calculate the imaginary component of 
the dielectric function and if desired the real component, 
using the appropriate Kramers-Kronig relation. In prac- 
tice, the joint DOS term is smeared using a Gaussian 
function, and the scissor operator is often applied to fa- 
cilitate comparison with experiment. 

One could envisage constructing polymer/Cgo devices 
in a variety of ways with either covalent or non-covalent 
bonds connecting the polymer to the Cgo, or by creat- 
ing thin film devices. Indeed there are a wide variety 
of structures which have already been realized for Cerj- 
polymer hybrids^. We have chosen simple structures for 
hybrids of Ceo with PPP and PPV with concentrations 
of approximately 25% by weight, in agreement with that 
used by Lane et al^ in their experiment on Cgo-dopcd 
PPP. We have calculated the optical absorption spectra 
for these hybrid structures, and compared the results to 
those of isolated chains of PPP and PPV and an isolated 
molecule of Cgo- 

Before generating structures for the hybrid polymers, 
individual calculations were performed for the two poly- 
mers and the fullerene, for validation. All calculations 
have been performed using the LDA exchange-correlation 
functional, with V point sampling only, norm-conserving 
pseudopotentials and periodic boundary conditions. The 
total energies were converged with respect to the psinc 
cut-off energy and NGWF radii - for all calculations the 
psinc grid spacing was set to be equivalent to a kinetic en- 
ergy cut-off of 1115 eV^i, valence NGWF radii of 5.29 A 
and conduction NGWF radii of 7.41 A were used for all 
atoms, with four valence and conduction NGWFs per C 
atom and one valence and four conduction NGWFs per 
H atom. No truncation was applied to the density kernel. 



III. RESULTS AND DISCUSSION 

In order to test the accuracy of the method compared 
to traditional plane-wave DFT, we have calculated the 
density of states and the imaginary component of the di- 
electric function for a small unit cell of PPV containing 
56 atoms, using ONETEP both with and without the op- 
timization of a set of conduction NGWFs. These results 
arc compared to PWPP results calculated using CASTEP 
in Figs. [5] and [3] respectively. The same pseudopoten- 
tials were used in both codes and the kinetic energy cut- 
off of the psinc and plane-wave basis sets were chosen 
to be approximately equivalent. As can be clearly seen 
from these results, the optimization of a set of conduc- 
tion NGWFs is vital for good agreement with traditional 
DFT results. Further examples demonstrating the accu- 
racy of the method are given elsewhere, together with a 
discussion on some of the limitations of the metho d 21 ' 22 . 

Two different bonding structures were used for Cerj- 
PPP, which are depicted in Figs. |4(a)| and |4(b)[ The C 60 - 
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FIG. 2. (color online) Density of states of PPV calculated 
with both ONETEP and a PWPP code. In the upper panel, 
the PWPP result is plotted with the ONETEP result in the 
valence NGWF basis only, in the lower panel the PWPP re- 
sult is plotted with the ONETEP result in the joint valence- 
conduction NGWF basis. For each result, a Gaussian smear- 
ing width of 0.1 eV is used and only fifteen conduction states 
are included, which corresponds to the number which were op- 
timized for the ONETEP conduction calculation. Reproduced 
with permission from Ref.—. 
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FIG. 3. (color online) The imaginary component of the dielec- 
tric function of PPV calculated using the PWPP method and 
ONETEP both with (lower panel) and without (upper panel) a 
conduction calculation. A Gaussian smearing width of 0.1 eV 
is used. Reproduced with permission from Ref.—. 



PPV was constructed using the same attachment mech- 
anism as structure A. Unit cells of 676 (672) atoms were 
constructed in periodic boundary conditions for the hy- 
brid polymers in structure A (B), with two molecules of 
Ceo evenly spaced across the polymers. Ground state and 
conduction calculations were then performed in ONETEP, 
followed by the calculation of the imaginary component 
of the dielectric function. The calculations on pure PPP 
and PPV were calculated in unit cells containing 280 and 
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(a)Stmcture A 




(b) Structure B 




(c)Structure B relaxed 

FIG. 4. Schematics showing segments of the two different 
structures of C60-PPP, including both the unrelaxed and re- 
laxed forms of structure B. 



560 atoms respectively, equivalent to samplings of 14 and 
20 /c-points with respect to the primitive unit cell, such 
that the spectra were converged. 

The structures of both PPP and PPV were optimized 
using CASTEP and the structure of Cgo was relaxed in the 
same manner. To reduce the computational cost, the hy- 
brid structures were initially not optimized. The forces 
on the structures were calculated in ONETEp2^; the root 
mean squared (rms) forces for structures A and B of the 
hybrid PPP were 0.4 and 0.5 eV/A respectively and the 
maximum forces were and 4.1 and 4.7 eV/A respectively. 
For the hybrid PPV they were 0.4 and 5.1 eV/A. For the 
separated structures the forces were less than 0.1 eV/A. 
As these forces for the hybrid structures are rather high, 
a test was done for structure B of Cgo-PPP, by optimis- 
ing the structure using ONETEP to give a rms force of 
0.02 eV, and a maximum force of 0.1 eV. The changes 
in the structure are essentially localized to the section of 
the chain closest to the fullerene, with the polymer be- 
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ing bent away from the Ceo, as shown in Fig. 4(c 



resulting spectrum was then compared to that of the un- 
optimized structure, giving an estimate of the accuracy 
of the results for the other structures (see Fig. [Sj . 

Results for the different hybrid PPP structures are 
plotted in Fig. [5J with the spectra of the PPP and Ceo 
shown for comparison. The sum of the separate poly- 
mer and Ceo results is also plotted. The results for both 
structures of Cgo-PPP show a number of small additional 
peaks below the original onset of absorption of PPP (ap- 
proximately 2.4 eV), which are not present in either PPP 



or Ceo on their own. These peaks are stronger for struc- 
ture B than structure A. The relaxation of structure B 
has had a noticeable impact on the spectra, with changes 
both in position and relative height of the observed peaks. 
In particular, those peaks which are below the original 
onset of absorption have significantly decreased in height, 
which is likely due to the bending of the polymer away 
from the fullerene. However, the overall picture of addi- 
tional peaks below this onset is still clearly preserved. 

The imaginary component of the dielectric function has 
also been calculated for Cgo and PPP in the same unit 
cell with the same cell dimensions as for the hybrid struc- 
tures, but separated by a large distance (greater than 
14 A) so that they are not chemically bonded, the result 
for which is also shown in Fig. [SJ In this case, the relative 
heights of the peaks have been affected, and an additional 
peak is present at around 3.5 eV, however no additional 
peaks are present below the original onset of absorption. 
This confirms that for those peaks below 2.4 eV to occur, 
the Ceo and PPP must be in close proximity. 

One of the advantages of computational spectroscopy 
is that it can be used to determine which transitions are 
responsible for a given peak. In this way insight can 
be gained into the origins of the additional peaks in the 
hybrid PPP system. For example, for the unoptimized 
structure B we can identify that there is a transition at 
2.05 eV between the HOMO-1 (highest occupied molec- 
ular orbital) and the LUMO+12 (lowest occupied molec- 
ular orbital), for which the wavefunctions arc plotted in 
Figs. 7(a) and 7(b) We can see that the electron den- 



sity is primarily located on the polymer for the HOMO- 
1, whereas for the LUMO+12 it is primarily located on 
the fullerene. However, the three additional peaks below 
2.4 eV result from a large number of transitions that are 
close together in energy. Therefore, rather than attempt- 
ing to visualize all of the individual transitions involved, 
as we have done for this example, it is more useful to gen- 
erate an average density change, weighted by the strength 
of the transition involved. To this end, we have used the 
following expression to calculate a weighted density dif- 
ference n we ightod: 



E v JMcv\ 2 (NAcI 2 



^-weighted 



E„,c We 



(5) 



where the sum is over all valence and conduction states 
which correspond to energy transitions in some given 
energy window and M cv is the optical matrix element 
for the transition between the corresponding valence and 
conduction state. Fig. 7(c) shows the resulting density 



for transitions below 2.4 eV for the unoptimized struc- 
ture B. As can be seen, there is a significant increase in 
electron density on the Cgo, whereas the electron density 
has decreased along the length of the chain. This agrees 
with the expectation from experiment that charge trans- 
fer will occur between PPP and Cgo and suggests that the 
charge transfer process is stronger for structure B of the 
hybrid polymer. This is possibly due to the fact that the 
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FIG. 5. (color online) The imaginary component of the dielectric function of C60-PPP in structures A and B, with both the 
unoptimized and optimized structure of B, compared to the sum of the spectra of pure PPP and Ceo- The spectra of pure PPP 
and Ceo are also compared with the total, as is that of Ceo and PPP in the same cell separated by a large distance (denoted Ceo 
with PPP in the key). The three insets show a closer view of the spectra between 0.5 and 2.5 eV, highlighting the additional 
peaks. The results are plotted with a Gaussian smearing width of 0.05 eV. 
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FIG. 6. (color online) The imaginary component of the dielec- 
tric function of the Ceo-PPV hybrid, compared to the sum of 
the spectra of pure PPV and Ceo- The spectra of pure PPV 
and Ceo are also compared with the hybrid spectrum in the 
lower half of the plot. The results are plotted with a Gaussian 
smearing width of 0.05 eV. 



Ceo is more strongly bonded to the chain in structure B, 
with six C-C bonds rather than two. More importantly, 
we can see both from the bending of the chain during the 
relaxation of structure B and by visual inspection of the 
cigenstates close to the band gap that the pi-bonds on 
the phenyl ring adjacent to the fullerene have been dis- 
rupted, allowing for a strong mixing between states - an 

For structure A, 



example of this can be seen in Fig. 7(a) 
on the other hand, the angle between the fullerene and 
the polymer is such that no similar disruption occurs, as 
shown for the HOMO-4 in Fig. |7(d)| As a result, the ad- 
ditional peaks in the absorption spectrum for structure 



A are much less pronounced than those of structure B. 

The spectrum for the hybrid Cgrj-PPV is shown in 
Fig. |Hl again plotted with those of PPV and Ceo- Unlike 
the PPP hybrids, the spectrum does not show any sig- 
nificant additional peaks, although some peaks have in- 
creased in intensity. The lack of additional peaks in the 
spectrum for Ceo-PPV could indicate that DFT alone is 
insufficient for capturing such behaviour in this particu- 
lar system - it has been previously demonstrated that the 
inclusion of excitonic effects has a strong impact on the 
absorption spectrum of PPV—. However, an alternate 
study by Pedersen using density-functional-bascd tight- 
binding (DFTB) indicates that excitonic effects are im- 
portant for both PPP and PPV— and so one also might 
not expect to see such peaks in Cgn-PPP without their 
inclusion, which we can see is not the case. Alternatively, 
it could be that the structure considered here is not ideal 
for the transfer of charge; in a further study by Peder- 
sen^ charge transfer states are calculated in Ceo-PPV 
with an arrangement similar to structure B, but with a 
finite chain length, a range of separation distances be- 
tween the Ceo and the PPV greater than that used here 
and again using DFTB with the inclusion of excitonic ef- 
fects. Given that both the structure and formalism are 
different from those used in our calculation it is difficult 
to say which might have the bigger effect, however as the 
low energy peaks in structure A of Ceo-PPP are already 
very small it seems more likely that it is the structure 
that is unfavourable, with too large an angle between 
the fullerene and the chain. We can also see from the 
spectra that the states on Ceo are separated in energy 
from those of PPV, whereas with PPP there is a signifi- 
cant overlap in energy, particularly around 2.6 eV. This 
should make it easy for these states to mix in the hybrid 
PPP structure to create a lower-lying state, which would 
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(a)HOMO-l structure B 




(b)LUMO+12 structure B 
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(c)Weightcd density difference for structure B 
(d)HOMO-4 structure A 

FIG. 7. (color online) Two different eigenstates for a segment of unrelaxed structure B of Ceo-PPP, one occupied and one 
unoccupied, a weighted density difference for transitions below 2.4 eV for structure B and an eigenstate of structure A. For the 
weighted density, red indicates an increase in electron density and blue indicates a decrease. These figures were plotted using 
XCrySDen^. 



not be the case for PPV. The creation of such a state 
seems to be confirmed by the significant increase in the 
absorption spectra just below this energy seen for both 
PPP hybrids and for PPP and Ceo in the same simula- 
tion cell. Finally, it could also be possible that there is 
some transfer of charge occurring between the PPV and 
the Ceo, but it is not visible in the absorption spectrum 
cither because it is too weak, or because it is masked by 
other peaks which are already present. 

In practice, it is likely that chemically substituted 
PPP and PPV or other related conjugated poly- 
mers, such as poly[2-methoxy, 5-(2-cthylhcxoxy)-l,4- 
phenylene vinylene] (MEH-PPV) or poly[2-methoxy-5- 
(3',7'-dimethyloctyloxy)-p-phenylene vinylene] (MDMO- 
PP V) , might be used instead of pure PPP and PPV, due 
to the increased solubility that can be achieved in organic 
solvents, which is useful for easy fabrication of electronic 
devices. Similarly, a lot of interest has been generated 
by fullerene derivatives such as [6,6]-phenyl-C6i butyric 
acid methyl ester (PCBM) as an alternative to pure Ceo 
- there have, for example, been a number of studies on 
MDMO-PPV in conjunction with PCBM, including the 
development of new techniques for experimental charac- 
terization^—. Furthermore, it would not necessarily be 
possible to easily manufacture the structures investigated 
in this work. The above simplifications are justified as 
a first approximation in allowing us to model the ab- 
sorption spectra of fullerene-conjugated polymer hybrids 



whilst keeping the computational costs low. However, it 
would be interesting to extend this work in future to the 
consideration of other hybrids such as those mentioned 
above, which, despite the larger unit cells involved, will 
still be fully accessible within a LS-DFT framework. 



IV. CONCLUSIONS 

In conclusion, we have investigated the effect on the op- 
tical absorption spectra of the attachment of molecules 
of Cgo to the two conjugated polymers PPP and PPV 
using a LS-DFT formalism, which has allowed us to sim- 
ulate long chains with a realistic proportion of Cgo with 
respect to the polymer chain length. For PPP, we con- 
sidered two different bonding structures between the Cgo 
and the chain. For both structures of the Ceo-PPP hybrid 
additional peaks were present in the absorption spectra 
below the original onset of absorption. By identifying the 
transitions between Kohn-Sham eigenstates which were 
responsible for these peaks and using the optical matrix 
elements to form a weighted density difference, we were 
able to demonstrate the transfer of electrons from the 
length of the chain to the Cgo molecule, agreeing with ex- 
perimental observations of charge transfer. For the PPV 
we did not find any additional peaks below the original 
onset of absorption, which is likely an indication that the 
bonding structure selected here is not optimal. Further 
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calculations could be done in future on different bonding 
structures of C 60 -PPV such as structure B, in order to 
confirm this theory. 

A number of approximations and simplifications were 
made, including the use of a purely DFT based formal- 
ism, which was justified based on the impracticality of 
going beyond DFT for such large systems. In addition, 
relaxation of the hybrid structure was shown to affect the 
details of the absorption spectrum but not the overall pic- 
ture and thus the use of unoptimized structures was also 
justified. Furthermore, whilst the structures used were 
not necessarily those directly achievable in experiment 
or likely to be the precise systems of interest for photo- 
voltaic systems given that interest has shifted from pure 
PPP and PPV to various derivatives, this work provides 
an ideal starting point for the study of large fullerene- 



polymcr complexes. Finally, this work could be easily 
extended in future to consider more complex polymers 
and fullerene-polymcr combinations; indeed the linear- 
scaling formalism used here is ideal for the computational 
study of such materials and should prove to be useful in 
gaining further insights into systems which are of great 
interest for organic photovoltaics. 
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